Magnetostructural Effect in the Multiferroic BiFe03-BiMn03 Checkerboard from 

First Principles 
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Using first principles calculations, we present a magnetostructural effect in the BiFeC>3-BiMn03 
nanocheckerboard that is not found in either bulk parent compound or in BiFe03-BiMn03 super- 
lattices. We also demonstrate that the atomic-scale checkerboard has a multiferroic ground state 
with the desired properties of each constituent material: polar and ferrimagnetic due to BiFeOs and 
BiMnC>3 respectively. 

PACS numbers: 75.80.+q,77.80.-e,75.75.+a 



There is currently tremendous interest in finding new 
multiferroic (ferroelectric and ferromagnetic) materials 
with large magnetoelectric coupling. Advances in the 
synthesis of artificially structured materials have stimu- 
lated efforts to design new multiferroic heterostructures, 
with first principles methods being an essential tool for 
the identification and investigation of promising systems. 
In this Letter, we report the first-principles identifica- 
tion and characterization of an unusual heterostructurc, 
a multiferroic atomic-scale 2D nanocheckcrboarcU^^ of 
BiFe03-BiMn03, with properties that critically depend 
on the geometry and are not present in either bulk or 
layered structures of the constituent materials. In par- 
ticular, the 2D checkerboard geometry leads to magnetic 
frustration and to quasi-degenerate magnetic states that 
can be tuned by an external perturbation that changes 
the crystal structure, such as an electric field. This re- 
sults in a novel magnetostructural effect, adding to pre- 
vious examples of magnetostructural coupling such as 
bulk 5 - and layered^ manganites, epitaxial EuTiOjL and 
EuSej ' PbSe\- x Te x multilayers^ 

Our first principles calculations are performed using 
density functional theory within the local spin-density 
approximation (LSDA)+U method as implemented in 
the Vienna Ab-initio Simulation Package VASP-4.6.342.. 
We test the robustness of our results with two different 
implementations of the rotationally invariant LSDA+U 
version, the first as introduced by Liechtenstein^ with 
Ufc = U]\i n = 5eV, Jfc = Jhin = leV, and the sec- 
ond due to Dudare^ii, with = 5.2eV", U e JJ = 4eV, 
where U e ^ — U — J. It has been shown that these U 
and J values match experimental data in bulk BiFeO^ 1 ^; 
the value U e ^ = 5.2eV has been used for previous bulk 
BiMnC>3 ground state calculations 1 ^. We use projector- 
augmented wave potentials (PAW) 14 ' 15 and treat ex- 
plicitly 15 valence electrons for Bi (5d 10 6s 2 6p 3 ), 14 for 
Fe (3p 6 3d 6 4s 2 ), 13 for Mn (3p 6 3d 5 4s 2 ), and 6 for O 
(2s 2 2p 4 ). The cutoff energies for the plane wave ba- 
sis set are 550eyand 800eV in the ionic relaxations and 
for subsequent self-consistent energy calculations respec- 
tively. Gaussian broadening of the partial occupancies 
for each wavefunction is 0.05eV. A Monkhorst-Pack 
k-point gricU^ is generated with density 4 x 4 x 4 for 
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Figure 1: (i) BiFe03-BiMn03 superlattice with alternation 
of Fe/Mn planes, (ii) (left) BiFe0 3 -BiMn0 3 checkerboard. 
Checkerboard ordering of Fe/Mn atoms in the (xy) plane, 
pillars of the same composition form along the z-direction. 
(right) Ideal perovskite unit cell. Perovskite cells with Fe/Mn 
atoms on the B-site repeat according to the checkerboard pat- 
tern (ii), or layered geometry (i). 



x s/2 x l)ao double perovskite and 4 x 4 x 2 for 
(\/2x y/2x 2)ao four perovskite cells. Ions are relaxed to- 
wards equilibrium positions until the Hellmann-Fcynman 
forces are less than 10~ 3 eV/ A. The spontaneous polar- 
ization is calculated by the Berry phase method^ with 
k-point mesh twice as dense as in the energy calculations. 

We consider four formula units (perovskite cells), two 
each with Fe and Mn atoms on the i?-site, which we re- 
peat periodically in space. For the planar checkerboard, 
we alternate iron (Fe) and manganese (Mn) atoms at the 
atomic level to form pillars of the same composition as 
in Fig. [IJii). For the layered superlattice, we alternate 
single unit cell layers along z, as in Fig. [Hi). In both 
cases, the supercell is V2clq x \/2ao x 2a,Q. 

We study various collinear spin orderings of the mag- 
netic Fe and Mn atoms, shown for the checkerboard in 
Fig. [21 FcFM and FeAFM refer to ferromagnetic and 
antifcrromagnetic ordering respectively for the Fe mo- 
ments in the relevant structural component (pillar for 
the checkerboard, layer for the superlattice); similarly 
MnFM and MnAFM describe the spin ordering of the 
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Figure 2: From top left to bottom right: (i) G-AFM: rocksalt 
type antiferromagnetic (AFM) order, (ii) C-FIM: AFM order 
in horizontal planes, ferromagnetic (FM) order along Fe/Mn 
pillars, (iii) FM order, (iv) FeAFMMnFM: AFM order along 
Fe pillars, FM order along Mn pillars, (v) FeFMMnAFM: 
FM order along Fe pillars, AFM order along Mn pillars, (vi) 
FeAFMMnAFM: AFM order along Fe/Mn pillars, but FM 
order in horizontal planes. 



Mn moments. In the checkerboard, this notation fully 
specifies the states considered. For the superlattice, the 
remaining ambiguity is resolved as follows: FeAFMM- 
nAFM magnetic order designates AFM ordered Fe and 
Mn planes with FM order along the mixed Fe-Mn chains 
in the z direction, while G-AFM designates the case with 
AFM order along the mixed chains; similarly, FcFMM- 
nFM designates FM ordered Fe and Mn planes with AFM 
order along the mixed Fe-Mn chains in the z direction, 
while FM designates the case with FM order along the 
mixed chains. 

In searching for the ground state crystal structure for 
each spin ordering, we consider structures generated by 
three typically unstable modes of the cubic perovskitc 
structured the zone center polar mode TJ, the M^" oxy- 
gen octahedron rotations (all rotations about a given axis 
in phase) and R|" rotations (sense of rotations alternates 
along the rotation axis). We freeze in selected modes 
and combinations of modes and optimize atomic dis- 
placements and lattice parameters in the resulting space 
groups. 

First, to investigate the effect of -B-site cation geometry 
on the magnetic ordering, we present results in Table U 
for the magnetic ordering energies when the structures 
are held fixed to the ideal perovskitc reference structure. 
In the layered superlattice and both bulk parent sys- 
tems, the difference in energy between magnetic ground 
state (FeAFMMnFM in the superlattice, G-AFM in bulk 
BiFe0 3 and FM in bulk BiMn0 3 ) and the first alter- 
native state is 0.11 — 0.12eV/ f.u.; this difference corre- 
sponds to a relatively large energy and we do not ex- 
pect a transition to a different magnetic state. The high- 
est energy magnetic states are more than 0.26eV/ f.u. 
apart. On the other hand in the checkerboard, all mag- 
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Figure 3: Structural energetics of bulk BiFeOg. Energy differ- 
ence per perovskite cell (f.u.) for different magnetic orderings 
(see Fig. [2} and structural distortions: (1) Pm3m: no dis- 
tortion - ideal perovskite, (2) Fj(z): polar distortion along 
z axis, (3) M^(z): + oxygen octahedra tilts about z axis, 
(4) Rf(y): — oxygen octahedra tilts about y axis, (5) Rfiy) 
and T^(y) (Rr(y)): linear combination of (4) and (2) along 
y axis, relaxes back to polar TJiy) with zero tilting angle, 
(6) .Rj ([111]) (Rt(d)): - oxygen octahedra tilts about [111] 
axis, (7) Rj([lll]) and FJ([111]) (RV(d)): linear combina- 
tion of (6) and (2) along [111] (d), where d refers to the cube 
diagonal direction. 



netic states are found quasi-degenerate and are confined 
within the energetical window of 0.12eV/ f.u., that is, all 
are lower than the lowest states in the layered super- 
lattice and bulk parent compounds. Indeed, the closest 
magnetic state to the FeAFMMnFM ground state is now 
only 0.022eV/ f.u. higher, making it much more plausible 
for a magnetic transition to occur. 

Next we study the energetics of the structural distor- 
tion and its effect on the spin order. Before discussing 
results for the BiFe03-BiMn03 checkerboard, we look at 
the structural energetics of the two bulk constituent ma- 
terials, BiFe03 and BiMn03. We plot energies for various 
magnetic orderings in seven types of structural distor- 
tions of bulk BiFcOs in Fig. [3] Our calculation verifies 
the R3c ground state of BiFeOs: counter-rotations of the 
oxygen octahedra (i?4~) and polar ionic distortions (rj) 
along the [111] axis are most energetically favorablei 12 i 19 
The ground state structure has G-type AFM order and 
spontaneous polarization 90/iC/cm 2 along [111] axis. For 
all structural distortions considered, the lowest energy 
magnetic ordering is G-AFM. 

We study the structural energetics of bulk BiMnO.3 
in a similar way; the plotted energies for various struc- 
tural distortions and magnetic orderings are presented in 
Fig. 21 We find the lowest energy structure with R3c sym- 
metry, the same as the ground state of BiFeOs. The low- 
est energy structure has a half-metallic character and is 
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Table I: Calculated magnetic energies in an ideal perovskite setting with lattice constant ao — 3.893A for various magnetic 
states in the checkerboard, layered superlattice, and bulk BiFeOs and BiMnOs. Value of U = 5eV and J = leV is used. 
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Figure 4: Structural energetics of bulk BiMnOs. Energy dif- 
ference per perovskite cell (f.u.) for various structural dis- 
tortions (see Fig. |3J) and magnetic orderings (see Fig. [2); 
P4/mmm corresponds to a tetragonally distorted perovskite 
cell with ideally positioned atoms and MF(z) is a linear com- 
bination of rotational M^(z) and polar T±(z) modes. 



Figure 5: Structural energetics of BiFe03-BiMn03 checker- 
board. Energy difference per perovskite cell (f.u.) for different 
magnetic orderings (see Fig.0 and structural distortions: (1) 
PA/mmm, (2) r7J», (3) R+(y) and F+(y), (4) Rt([lU]) and 
rj([lll]) (R± , rj (d)). Inset: zoomed view of the magnetic 
energies of c-R3c (4) distortion. 



ferromagnetic. This structure is not the monoclinic cen- 
trosymmetric ground state C2/c of bulk BiMnOs which 
has a larger unit cell than that considered here.— How- 
ever our calculation shows that it lies close to the ground 
state (only 43 meV/f.u. above the GS). For all struc- 
tural distortions considered, the lowest energy magnetic 
ordering is FM. 

In the layered BiFc03-BiMn03 superlattice, we cal- 
culate magnetic energies for the rocksalt type G-AFM 
and FcAFMMnFM layered magnetic states in two struc- 
tural distortions. For (y)&,Y^ (y), we find AE = 
-0.504eV/f.u. for G-AFM and AE = -0.553eV/ f.u. for 
FeAFMMnFM with respect to the FeAFMMnFM mag- 
netic state in the ideal perovskite cell (see Table H]). For 
i? 4 f ([lll])&r 4 -([lll]), we find AE = -0.752eV/ f.u. for 
G-AFM and AE = -0.76W/ f.u. for FeAFMMnFM. 
For both structural distortions considered, the lowest en- 
ergy magnetic ordering is FcAFMMnFM. 



Let us now look at the results for the structural en- 
ergetics of the BiFc03-BiMn03 checkerboard. In Fig. 
we present the energies for four different types of struc- 
tural distortions. These types of distortions show the 
lowest energies among a larger set of structures that we 
explored3i Notice that the variation of the structural en- 
ergy is much larger than that of the magnetic energy of 
the checkerboard. 

Not surprisingly, the P+([lll]) and 1^" ([111]) (R3c) 
type of structural distortion is energetically the most fa- 
vorable; it is the BiFe03 ground state and the BiMnOs 
lowest energy structural distortion. The R3c symmetry 
is now broken due to pillar cation ordering and the space 
group of the checkerboard ground state becomes PI; 
we use the notation c-R3c, where c designates "checker- 
board", as a reminder of the origin of the distortions. 
As shown in the inset of Fig. the two lowest magnetic 
states G-AFM and FeAFMMnFM, arc only 2meV/f.u. 
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apart. The ground state of the checkerboard has the 
FeAFMMnFM magnetic order, where Fe spins are or- 
dered antifcrromagnctically (AFM) along the Fe pillars, 
Mn spins are ordered ferromagnetically (FM) along the 
Mn pillars, reflecting "AFM" and "FM" nature of the 
parent BiFeC>3 and BiMnC>3 compounds respectively. A 
total magnetic moment 3.7 [Ib per Fe -Mn pair results 
from manganese chains. The FeAFMMnFM ground state 
is insulating with energy gap 0.88eE, and we calculate 
a value of the polarization 62 u\C/ cm 2 pointing in the 
[0.85,0.85, 1] direction. The ground state of the checker- 
board is multiferroic, being ferroelectric and ferrimag- 
netic. 

In particular we want to relate and contrast the prop- 
erties of the BiFc03-BiMn03 checkerboard to those of its 
two bulk constituent materials; we recall that BiFe03 is 
polar and antiferromagnetic, while BiMnO"3 is non-polar 
and ferromagnetic, and we have found that the checker- 
board assumes the desired ferromagnetic-ferroelectric 
properties of each leading to a multiferroic ground-state 
whose magnetic behavior is structurally sensitive. We 
attribute this behavior to the development of a quasi- 
degenerate manifold of magnetic states in the checker- 
board, in contrast to the gapped states in the layered 
superlattice and in the bulk; this can be understood in 
terms of frustration of the cations inherent in the checker- 
board geometry. Since bulk BiFeC>3 is known to be G- 
AFM, and bulk BiMn0 3 FM, the Fe-Fe and Mn-Mn 
interactions tend to be AFM and FM-like respectively. 
In the layered superlattice, each Fc(Mn) atom has four 
Fe(Mn) and only two Mn(Fe) nearest neighbours, so that 
the Fe/Mn layers prefer to be AFM/FM, leading to mi- 
nor frustration between the minority of mixed Fe-Mn 
bonds. The FeAFMMnFM layered ground state is much 
more preferable and lower in energy than any other mag- 
netic state. In the checkerboard there are more frustrated 
bonds per each cation, and therefore more weight is given 
to the mixed Fe-Mn bonds and various magnetic states 
compete energetically. 



We study the sensitivity of the closely spaced magnetic 
levels in the checkerboard to a structural distortion. As 
we tune the system from the checkerboard c-R3c ground 
state to c-I4cm state with (y)kT^ (y) distortions, ei- 
ther the FeAFMMnFM (filled diamond) or the G-AFM 
(open circle) lowest state is favored. Switching between 
these two magnetic states occurs as we tune the sys- 
tem to other structural distortions (see Fig. . It is the 
competition between these two magnetic types that al- 
lows switching between nonzero and zero magnetization; 
the magnetostructural effect leads to the possibility of 
a structurally-driven magnetic transition in the checker- 
board. This could be realized, for example, by imposing 
cxpitaxial strain constraints. 

In summary, we present a magnetostructural effect in 
the atomic-scale checkerboard BiFe03-BiMn03, which 
is not present in either bulk or in layered structures of 
these two materials. Furthermore, unlike its parent com- 
pounds, the checkerboard has a multiferroic ground state 
with a nonzero magnetization and polarization; this is a 
new example of a nanocomposite with properties that 
can be designed. We note that this behavior is due to 
the magnetic frustration in this system inherent to the 
checkerboard geometry; as a result the magnetic states 
are quasi-degenerate and can be tuned by small perturba- 
tions including strain. We remark that our first principles 
calculations do not include spin-orbit coupling which is 
known to lead to weak ferromagnetism in bulk BiFe03^ 
We expect that such corrections will not change our re- 
sults fundamentally, but this is certainly worth pursuing 
in future work. We would also plan to investigate simi- 
lar checkerboards on longer length-scales to make more 
direct contact with the possibility of future experiments. 
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